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Abstract. We extend the results [19], [18], [17] about the fluctuations of the 
matrix entries of regular functions of Wigner matrices to the case of sample 
covariance random matrices. 



1. Introduction and Main Results 

Recently, there have been a number ol results concerning matrix entries of func- 
tions of random matrices. That is, for a N x N random real symmetric (Hermitian) 
matrix, Mjsr, we consider the entries of the matrix /(Af/v) where / is a regular test 
function. 

In [14], Lytova and Pastur consider the case where Mpj is drawn from the Gauss- 
ian Orthogonal Ensemble (GOE) or Gaussian Unitary Ensemble (GUE). We recall 
that a GOE matrix is defined as M N — -±=(Y N + Y N ), where the entries of Y N are 

i.i.d. iV(0, \cr 2 ) real random variables (see e.g. [1]). In a similar way, a GUE matrix 
is defined as Mm — — L(Yftr + Y N ), where the entries of Xjy are i.i.d. N(0, jcr 2 ) 
complex random variables. It was shown in [14] that 

Vn(J{Mn)h -E[/(M w ) y ]) — ► N ( ,i±^ 2 (/)) , (1.1) 

in the limit when the size of the matrix goes to infinity, where w 2 (/) = Y(f(ip)), ip 
is a random variable distrubited according to Wigner semicircle law, and /3 = 1 for 
the GOE and j3 = 2 for the GUE. We recall that the Wigner semicircle distribution 
is supported on the interval [— 2cr, 2cr] and its density with respect to the Lebesgue 
measure is given by 

^(x) = -l^y/A^-x^^x). (1.2) 

In the case where Wm = ~j n °An, arid is a symmetric (Hermitian) Wigner 
matrix ([1], [5]) with i.i.d. (not necessarily Gaussian) entries up from the diagonal, 
Pizzo, Renfrew, and Soshnikov studied in [19] the fluctuations of both the diagonal 
and off-diagonal entries under the condition that the off-diagonal entries of An 
are centered and have finite fourth moment, and the diagonal entries of An are 
centered and have finite second moment. The variance of the off-diagonal entries, 
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as before, is equal to a 2 . The test function / has been assumed to be four times 
continuously differentiablc. In particular, it is shown in [19] that 

(f(W N )ij - E [f(W N )ij]) 

converges in distribution to the sum of two independent random variables: the first 
(up to scaling) is given by {A^)ij and the second is a Gaussian random variable with 
mean zero and variance explicitely given in terms of the function /. In addition, 
it was proven in [19] that the joint distribution of any finite number of normalized 
matrix entries converges to the product of one-dimensional limiting distributions. 
If the marginal distribution of the entries of Wn is Gaussian (so Wn belongs to 
the GOE (GUE) ensemble), one recovers (1.1). 

Such results might be considered as an analogue of the E.Borel theorem for 
the matrix entries of random matrices from the classical compact groups (see e.g. 
[8], [11], and [12]). In addition, the results about the fluctuation of the resolvent 
quadratic form are related to the limiting distribution of the outliers in the spectrum 
of finite rank deformations of Wigner matrices (see e.g. [20] and references therein). 

Almost simultaneously with [19] and using a different set of ideas, Pastur and Ly- 
tova [18] gave another proof of the limiting distribution of the normalized diagonal 
entries _ 

VN(f(W N )u-E[f(W N )u]), l<i<N 
when Wn — ^/n^-n, and An is a real symmetric Wigner matrix with i.i.d. cen- 
tered entries up from the diagonal provided the cumulant generating function 
log(Eexp(zA 12 )) is entire (so, in particular, all moments of the marginal distri- 
bution are finite) and the test function / satisfies 

/ (l + 2|fc|) 3 |/(fc)|dfc < oo 

JR 

where / is the Fourier transform 

/» = 4= / e- ikx f(x)dx. (1.3) 

The results of [19] and [18] are extended in [17] to the case of a Wigner matrix 
with non-i.i.d. entries where it was assumed that the off-diagonal entries have uni- 
formly bounded fourth moments, diagonal entries have uniformly bounded second 
moments, and certain Lindeberg type conditions for the fourth moments of the 
off-diagonal entries and the second moments of the diagonal entries are satisfied. 
The test function f(x) is assumed to satisfy 

/ (l + 2\k\) 2s \f(k)\ 2 dk<^, 

JR 

for some s > 3. 

In this paper, we study the fluctuations of matrix entries of a sample covariance 
random matrix. Namely, we consider the case where 

M N = j^A N A* Nl (1.4) 

and An is an N x n rectangular matrix with independent entries. We begin with 
some definitions. 

Definition 1.1. Let An = {(AN)ij) 1<i<N -i<j<n be an iVx n matrix with complex 
entries. We say the matrix An satisfies condition CI if 
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(i) {Re(^4jv)ij, Im(^4jv)ij : 1 < i < N; 1 < j < n} is a collection of independent 
random variables, 

(ii) each entry (Ajy)ij has mean and variance a 2 , 

(iii) each entry satisifies E(Ajv)? • = 0, 

(iv) sup^j-EK^jv)^-! 4 = m 4 < oo, 

(v) the entries satisfy the Lindeberg condition for the fourth moments, that is, 
for all e > 0, 

ij2n(ANh\ 4 i {l{AN)ijl>eVN} ^o (i.5) 



N 2 



as iV — y oo . 



Definition 1.2. Let A^v = (( j 4jv)ij)i<j<Ar-i<j<n be an JV x n matrix with real 
entries. We say the matrix An satisifies condition C2 if {(Ajv)ij : 1 < i < N; 1 < 
j < n} is a collection of independent real random variables and conditions (ii), (iv), 
and (v) hold from Definition 1.1. 

We define Xjy :— -±=An and Mn := XnX* n . Throughout this paper, we assume 
that cn ■— n/N — > c e (0, oo) as TV — > oo. 

Definition 1.3. Let BbeanNxN self-adjoint matrix with eigenvalues Ai, . . . , Ajv- 
The empirical spectral density of B is given by 

1 N 

i=l 

The limiting empirical spectral density of Mjy is known as the Marchenko-Pastur 
Law (see [3], [16]). 

Theorem 1.4 (Marchenko-Pastur). Suppose that for each N , the entries of An are 
independent complex (real) random variables with mean and variance a 2 . Assume 
n/N — > c e (0, oo) and for any e > 

^J2 E \( A "h\\\(A N)ij \>eVN}^0 (1-6) 



as iV — y oo . TTierc with probability one, the emperical density um n tends to the 
Marchenco-Pastur distribution, u a , c , with ratio index c and scale index a 2 where 



dx \ 0, otherwise, 

with a point mass at with weight (1 — c) when c < 1, and where 

u+ :=o- 2 (l + Vc) 2 , 
u_ :=a 2 (l-V~c) 2 . 

Remark 1.5. We note that the Lineberg condition (1.6) is implied by the Linde- 
berg condition for the fourth moments (1.5). 

Given a probability measure a on the real line, its Stieltjes transform is given by 

du(x) 



j 



z e C \ supp(^). 

jj z — x 
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For Im z ^ 0, we have the following bound for the Stieltjes transform of any prob- 
ability measure on M 



/ 

■n: 



d/i(x) 



z — X 



< 



1 



\lm{z)[ 



(1.7) 



The Stieltjes transform of /x CT)C is denoted by g a , c and is characterized as the solution 
of 

za 2 g^ c {z) + (a 2 (c - 1) - z)g a . c {z) + 1 = (1.8) 

that decays to zero as z — > oo. 

The Stieltjes transform of the expectation of the emperical spectral distribution 
of M N is given by 

gN{z) = E [ = E [tr N (R N (z))} 

JR Z — X 

where tr^v := -^Tr is the normalized trace and Rn{z) ■— (zL N — M N )~ X is the 
resolvent of M N . If it does not lead to ambiguity, we will use the shorthand notation 
Rijiz) for (R N (z))ij, l<i,j<N. 

For s > 0, we consider the space H s consisting of the functions : R — > IR that 
satisfy 



110112== (l + 2|fc|) 



2.s 



2 dk < oo. 



(1.9) 



We recall that C k (X) denotes the space of k times continuously differentiate func- 
tions onXcl and define the C k (X) norm 

d l f{x) 



max 



wnc(x) 

We now present our main results 



dx l 



, xex, o < z < fc . 



(1.10) 



Theorem 1.6. Let be a N x n random matrix with real entries that satisifies 
condition C2. Let m be a fixed positive integer and assume that for 1 < i < m 



I)l4 



(i) := lim — E\Ai 



exists and for all e > 
1 ™ 

— ^E|(AAr) 4i | 4 l { | (Ajv)i3 | >eAr i/4 } 



N-. 

3 = 1 

as N — > oo. Assume Cjv — > c € (0, oo) as N — > oo 
Then we have the following: 
(i) The normalized matrix entries 



0, l<i<m, (1.11) 
Zet f £ H s for some s > 3. 



{^(/(Mjvjy - E [/(Mat),,]) : 1 < i < j < m} 

are independent in the limit N — > oo. 
(wi) for 1 < i < j < m, 

y/N{f{M N ) ij -E[/(M JV ) y ]) —J- 7V(0, W 2 (/)) 
m distribution as N — > oo, wZiere 

w 2 (/) = v(/fac)). 



(1.12) 
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and rj c is a Marchenko-Pastur distributed random variable with ratio index c 
and scale index a 2 
(in) For 1 < i < m, 



<N(f{M N )u - 
in distribution as N 

and 



oo, where 



p(f) = E 



f(Vc) 



K 4 (i) := 777,4(1) — 3cr 4 . 



(1.13) 



(1.14) 



Theorem 1.7. Let An be a N xn random matrix with complex entries that satisi- 
fies condition CI. Let m be a fixed positive integer and assume that for 1 < i < m 

m 4 (i) := lim -VeU,,! 4 

j 

exists and for all e > (1.11) holds as N — > 00. Assume cn — > c e (0, 00) as 
TV — > 00 and let f G H s for some s > 3. TTien ?i>e ftcwe i/ie following: 
(i) The normalized matrix entries 

{ ViV {f(M N )ij - E \f(M N )ij]) : 1 < i < j < m} 

are independent in the limit N 
(ii) For 1 < i < j < m, 



00. 



N {f{M N )ij - E [f(M N ) 



N(0^ 2 (f)) 



in distribution as N — > 00 where N (0,oj 2 (/)) stands for the complex Gaussian 
random variable with i.i.d. real and imaginary parts with variance ^oj 2 (f) and 
u>(f) is defined in (1.12). 
(in) For 1 < i < m, 



N(f(M N ) u - E[f(M N ) u }) -+ N (o,uj 2 (f) + ^-p 2 (f)) 

in distribution as N — > 00 where to(f) is defined in (1.12), p(f) is defined in 
(1.13), and 

Ki(i) := m 4 (i) — 2er 4 . 

Remark 1.8. The limiting distribution of an entry in the sample covariance case 
is Gaussian and differs from the Wigner case ([18], [19]) where the limiting distribu- 
tion is given by a linear combination of an independent Gaussian random variable 
and the corresponding entry of the Wigner matrix. However, in the square case 
(c = 1) the limiting distribution of \/N (/(Mjv)ij — E [/(Mjv)^]) coincides with the 



limiting distribution of (ff(Wjv)ij — E [s^aO^'D > where g(x) = f(x ) and W\ 



N 



is a Wigner random matrix. This is not surprising since Mn is the NxN upper-left 
corner submatrix of the (N + n)x (N + n) matrix Z N n , where the NxN upper- left 
and n x n lower-right corner submatrices of Zn,u are both zero, the N x n upper- 
right corner submatrix of Zn,u is given by Xn, and the n x N lower-left corner 
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submatrix of Zn,u is given by X N . The limiting spectral distribution of Zn.u in the 
case n/N — > c = 1 is given by the Wigner semicircle law and the technique of [19], 
[17] in the square case can be extended without any difficulties to Zn.ji- 

_ 2 

Remark 1.9. The functions 1 and x ^° 2 are the first two orthonormal polynomials 
with respect to fi a ^ c (dx). Therefore, by the Bessel inequality, the variance of the 
limiting Gaussian distribution for the diagonal entries is zero if and only if the test 
function is linear and the marginal distribution is Bernoulli. For the off-diagonal 
entries, it immediately follows from (1.12) that the variance is zero iff the test 
function is constant on the support of the Marchenko-Pastur law. 

Remark 1.10. It follows from Proposition 4.1 and Lemma 2.2 that if / € C 7 (R) 
for the diagonal entries i = j (/ € C 6 (M) in the off-diagonal case i ^ j), one can 
replace E[f(M N )ij] in Theorems 1.6 and 1.7 by Sij J f{x)dfjb (TtCN {x). Moreover, as 
shown in Proposition 4.1, if supp(/) n R+ is compact, where R + = [0, oo) and / 
has seven continuous derivatives, then E[/(Mjv)»»] = / f(x)dfi a , CN (x) + (^) . If / 
has six bounded continuous derivatives on then E [/(Mjv)ij] = O (^) , i ^ j. 

We divide the proof of Theorems 1.6 and 1.7 into several sections. In Section 2, 
we apply a standard truncation lemma to the matrix entries of A n- Section 3 is 
devoted to computing the expectation and variance of the entries of the resolvent, 
Rn(z) 7 and Section 4 extends these results to more general functions. In Section 5, 
we prove a central limit theorem for entries of /(Mjv) where f(x) is a finite linear 
combination of the functions (z — x)^ 1 , zeC\M. Finally, we extend this result 
to more general test functions / € H s by an approximation argument. 

2. Truncation and Extremal Eigenvalues 
We note that by (1.5), we can choose a sequence — > such that 

^Y, E \( A »h\\\(A N)ij \>e N SN}^Q (2-1) 

N i,j 

as N — > oo. 

Lemma 2.1. Assume that An is an N x n matrix that satisifies condition CI in 
the complex case (condition C2 in the real case). Then there exists a random N xn 
matrix An with independent entries and a sequence ejv which tends to zero as N 
tends to infinity such that 

(i) the entries (Ax)ij have mean zero and variance a 2 , 

(ii) sup i:/ \(A N )ij\ < e N y/N, 
(Hi) sup Nl jEKA^iA 4 < oo, 

(iv) ¥{A N ^ A N ) — > as N -> oo. 

Proof. We present the proof in the case where the entries of An are real. The 
complex case follows a similar argument. We begin by selecting a sequence ejy — > 
such that (2.1) holds. Then let 

(A N )ij = (^Jv)ijl{| (Ajv ) y |< ejv VJv}- 

Define 

mNij — E( J 4jv)ij 
v 2 Nij =a 2 -E(A N f ij . 
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Then wc have that 



i \ <V \ (A \ n E|(A w ) I3 | 4 i {|(j4w)ij|>ewV1 y } / i 

\m Nij \ < \H\(A N ) l:j \l {KANhA>fNVW} < £ 3 ^3/2 - ° N 3/2 

N (2-2) 



and similarly 



V 2 Nij < E \( A N)ij\ 2l {\(A N ), j \>e N ^/N} < O (^) , (2.3) 

\m Nij \ < J^. (2.4) 

We now define (Ajv)ij to be a mixture of 

(1) (A N ) l3 with probability 1 - - and 

2 

(2) a Bernoulli random variable £ Ni j with probability 4. 

where we denote the mean and the second moment of £jvi? by (i Nij and rjy^ . 
We can choose £jvij such that 

(i) E(A N )ij = 0, 

(ii) VCijvJy = a 2 , 

(iii) (^4jv)ij < CcnVN for some absolute constant C. 

We now verify that such a construction is possible. Essentially, we have to show 
that one can choose ^Nij in such a way that (i) and (ii) are satisfied and 

\Hmj\ < C x e N ^N and T 2 Nij < C 2 e 2 N N (2.5) 

for some absolute constants Ci, C 2 > 0. Indeed, if this is the case, we can construct 
£,Nij — ^Nij +*pNij where ipNij is a symmetric Bernoulli random variable satisfying 
IVWijl < Ccat-v/TV where C is an absolute constant that depends on C\ and C2. 
This would immediately follow from (2.5). To verify (2.5), we note that 




Solving for ^Mij in the first equation and rjvij in the second and applying (2.2)-(2.4) 
yields the required bounds (2.5), verifying the claim. 

We note that without loss of generality we may assume C — 1 by our choice of 
the sequence ejv- 

Next by (2.2) and (2.3), we have that 

EKi^f < EK^vM 4 + + ^ (cnVF)* < m 4 + 0(1). 



8 



S. O'ROURKE, D. RENFREW, AND A. SOSHNIKOV 



To complete the proof of Lemma 2.1, we apply (2.2), (2.3), and (2.1) to obtain 



V(A N ^A N )< 



( \m Ni j\ 



UnVn 



+ 



J Nij 



-JV 



N 



< 



£%n 2 ^ 



E\{A N )i 



+ V(\(A N ) lJ \>e N VN) 



{|(-4iv)i 3 l>eNVlV} 



as N — »■ oo. 



□ 



We can now apply Lemma 2.1 to obtain a result on the norm of the matrix 
jfA N A* N . This result follows from [3, Theorem 5.9]. We present it here for com- 
pleteness. 

Lemma 2.2. Under the assumptions of Lemma 2.1, we have that \\-^A N A* N \\ — > 
a 2 (l + Vc) 2 * n probability as N — > oo. 

Proof. Since Theorem 5.9 from [3] does not apply directly to j^A^A* N , we simply 
note that by Lemma 2.1, it is enough to show || j^AmA* n \\ — > c 2 (l + \/c) 2 in 
probability. Theorem 5.9 from [3] now applies to the matrix -^AnA* n to obtain 



> (T 2 (l + y/c) 2 + X 







as N — > oo for all x > 0. The proof is then complete by noting that Theorem 1.4 
implies that, with probability 1, 



lim sup 

JV— >oo 



-^A N A* N 



>a 2 {l + sfcf. 



We also note that by (1.11), we can choose a sequence e^v 
1 



: ^ E l( AA, )^| 4l {l(^™), 3 |>e N JVi/4 } 



□ 

-> such that 

(2.6) 



as N — > oo for 1 < i < to. 

Lemma 2.3. Lei be a N x n complex (real) matrix that satisifies condition CI 
(C2) and (2.6) /or 1 < i < m, where to is a fixed positive integer. Then there exists 
a random Nxn matrix An with independent entries and a sequence ejy which tends 
to zero as N tends to infinity such that 
(i) (j4jv)ij = {Ajq)ij for to < i < N and 1 < j < n, 

(ii) the entries {A^)ij have mean zero and variance a 2 , 
(Hi) sup J - 1 < 1 < m \{A N )ij\ < ejyN 1 / 4 , 

(iv) supjv^ i:( -E|(Ajv)y| 4 < oo, 

(v) ¥(A N ^ A N ) — > as N -> oo, 

M H E,( E I^| 4 - E|i y | 4 ) -»• 0, 1 < i < to. 

The proof of Lemma 2.3 is very similar to the proof of Lemma 2.1 and the details 
are left to the reader. 

By Lemmas 2.1 and 2.3, we will assume that all the entries of An are bounded 
by 6nV~N and that the entries satisify conditions (ii)-(iv) of Lemma 2.3 for the 
remainder of the paper. Indeed, since the truncated matrix coincides with the 
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original with probability going to 1, it is enough for us to prove Theorems 1.6 and 
1.7 for the truncated matrix. 

We will also need the following lemma for controlling the expectation of the norm 
of X^X^j. 

Lemma 2.4. For any k > 1, there exists a constant C > (depending only on a, 
c, and k ) such that 

E[\\X N X* N \\ k ] < C 

for N sufficiently large. 



Proof. For any e > 0, 



E \\\X N X 



N 



E 



\XmX 



L {\\X N X* N \\<^{l+^+e} 



+ E 



\X N X%\ 



\\X N X* N \\>a2(l+VB) 2 +e} 



< [o- 2 (l + y/h) 2 + e] k + kj™t k - 1 F(\\X N X% 
By [3, Theorem 5.9], we have that 



> 



t + a 2 (l + VC) 2 ^jdt. 



jf 



d,t 



-fe-2 



t k ~ x ¥ (\\x N x* N \\ > t + <t 2 (i + Vc) 2 

<C'N- k - 2 £ t k ~ x (a 2 (l + VC) 2 +t-ej " ' dt = 0(^ fe - 2 ) 

for some constant C > 0. Thus, 

E[||X^|| fe ] < [a 2 (l + Vc) 2 +e] k + 1 
for N sufficiently large. □ 

3. Mathematical Expectation and Variance of Resolvent Entries 

This section is devoted to the estimates of the mathematical expectation and the 
variance of the resolvent entries. Throughout the section, we will consider the real 
case. The proofs in the complex case are very similar. It follows from Lemmas 2.1 
and 2.3 that for the purposes of the proof of Theorems 1.6 and 1.7 we can assume 
that An satisfies properties (i)-(iii) in Lemma 2.1 and properties (ii)-(iv) in Lemma 
2.3. Indeed, such a truncated matrix coincides with Am with probability going to 1, 
and, therefore, if the results of Theorems 1.6 and 1.7 hold for the truncated matrix, 
they also hold for Am- 

We begin by recalling the basic resolvent identity 

(zl - Ai)- 1 = (zl - A^- 1 - (zl - A 1 )-\A l - A 2 )(zl - A^ 1 (3.1) 

which holds for all z € C where (zl — A\) and (zl — A 2 ) are invertible. 

We will also use the decoupling formula (see for example [13] and [15]): for any 
real-valued random variable, £, with p + 2 finite moments and (j) a complex- valued 
function with p + 1 continuous and bounded derivatives the decoupling formula is 
given by: 

%(0) = E^% W (0) + E (3-2) 

where Hid &rc the cumulants of £ and e < Csup t |<^p +1 ) (t)|E(|£|P+ 2 ), C depends 
only on p. It follows from the proof of the decoupling formula in [15] that if |£| < K 
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with probability 1, then the supremum in the upper bound for the error term can 
be taken over t € [—K,K]. 

Recall that we denote the entries of the resolvent Rn{z) — {zIn ~ Mn)~ x of 
Mjv by Rij(z), 1 < i, j < N. Using (3.1), we can compute the derivatives of the 
resolvent with respect to any entry 

dR 

Q ' J = Rik{X* R)ij + (RX)uRkj- (3.3) 
Cam 

We now use (3.2) and (3.3) to compute the expectation and variance of the 
resolvent entries. 

Proposition 3.1. Let Mn = jjAjyA* N be a random real (complex) sample covari- 
ance matrix satisfying condition C2 (CI) and Rn{z) — {zIn — Mn)^ 1 ■ Then 

' P 6 {\lm{z)\-') 



E[R 4i (z)]=g (T , CN (z) + 0^ v . n 

for 1 < i < N, uniformly on bounded sets ofC\~. 



E[Rik(z)]=0 



P 5 (\Im(z)\-') 



N J (3.5) 

for 1 < i ^ k < N, uniformly onzeC\l 

, „ n ( p 4 (i Mgrog [p w (\\x N \\)(\\R N ( Z )r + ii R N (z)fi 2 )] \ 

>^=°{ N J (3.6) 

for 1 < i,j < N, uniformly onzeC\R. 

Here and throughout the paper Pk denotes a polynomial of degree k with non- 
negative coefficients. 

In (3.6) we have included the norm of the resolvent in the error estimate. This 
will be useful in the proof of Proposition 4.6. 

Proof. The following inequalities will be useful in our calculations: 

JV 

\Rik{z)\ <\Im(z)\-\ J2\ R ^)\ 2 ^ II^MH 2 < IM*)r 2 (3-7) 

3=1 

We first prove (3.4) and (3.5). We define the following sets on the complex plane. 
Let T be an arbitrary large number. Let L be a sufficiently large constant, to be 
chosen later. 

Q N := {z : \z\ < T + 1 and | lm(z)\ > LN~ 1/5 } O n := {z : | Im(z)| > £7V~ 1/4 } 

(3.8) 

Note that if z € Q C N n {\z\ < T + 1} then | Im(z)| 5 < L 5 iV _1 . When combined 
with (1.7) this implies 

nR lt (z)]-g a , CN (z)\ < 2|Im(z)r 1 =0^ |Im ^ )r ^ . (3.9) 

Similarly, if z e 0% then 

nRMW < \iHzT 1 = o ( |Im ^ )r5 ) • (3-io) 
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For the remainder of the proof of (3.4) and (3.5) we will assume that z is in Qj^ 
and Ojsf, respectively. 

The proof of both statements begins with the resolvent identity (3.1), and then 
an application of the decoupling formula (3.2). 

N n 

zE[R lk (z)} =S ik + J2Yl ®[Rij(z)XjlXkl] 

3=1 1 = 1 

2 2 

= S ik + a 2 E[R ik (z)tv N (R N (z)XX*)] + ^E[(R N (z)X X* R(z)) ik ] + -^-E[R ik (z)} + r N , 

(3.11) 

where rN is the third cumulant term coming from p — 2 and the error from trun- 
cating at p = 2. 

From the definition of the resolvent we have Rn(z)(zIn — XX*) = In, which 
implies Rn{z)XX* = zRn(z) — In- Applying this identity to (3.11) yields: 

zE[R ik (z)} =S ik + a 2 zE[R lk (z)tr N (R N (z))} - o- 2 E[R lk {z)] + 

a 2 no 2 ( 3 - 12 ) 

—E[(R N (z)XX*R N (z)) ik } + —E[R ik (z)\ + r N - 

We begin with the following lemma: 
Lemma 3.2. For z e C \ M: 

c^^MRMm < ( 3,3) 

r„ < (3.14) 
Additionally, for z £ On: 

rN < P 2 (|Im(z)|- 1 )E[P 8 (||A- jv ||)(||i 1 ' jv (z)|| 2 + \\R N {z)fl 2 ) ^ ^ 

Proof. To prove (3.13) we begin with the following bounds from Proposition 4 in 
[22]:' 

vrtr (r ( Z )))< |Im(zr4 vrtr m r~vn - I M*)r 7/2 E[l|iM*)ll] 3/2 

(3.16) 

It follows from the proof of Proposition 4 in [22] that these bounds are valid provided 
the fourth moments are uniformly bounded ([23]). Additionally, from (3.7) we have 

V(Rij{z)) < llm^r^Elp^z)!!] 3 / 2 . (3.17) 

Using Cauchy-Schwarz this implies 

as desired. 
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Now we prove (3.15); the argument along with Lemma 2.4 can be modified to 
prove (3.14). The third cumulant term in the decoupling formula is: 

j=i i=i L kl 

N n 

= 5Z^^((A A r) fc; )[i? 4 fc(2)(X*i?(z)X) H i 1 , fc ,(z)^ 
3=1 1=1 
1 ™ 

= Jjzji K 3 ((A N ) kl )E[R ik (z)(X*R N (z)X) ll (R N (z)X) kl 

+ {R N {z)X) u {R N (z)X) kl (R N (z)X) kl +R ik (z){R N {z)X) kl 

+ R ik (z){X*R N {z))i k (X*R N (z)X) H + {R N {z)X) a R kk {z){X*R N {z)X) u ], 

where ^((An)^) is the third cumulant of (A^ki- By condition C2, the «3((Ajv)fc/)'s 
are uniformly bounded. 

Using (3.7) and the Cauchy-Schwarz inequality this term is seen to be 

\im(z)\-'E[p 4 (\\x N \\)\\R N (z)\\ 2 } 



o 



N 



The truncation error is bounded from above by a finitely many sums of the 
following form 

^f]supE|i?: fc (z)(i?^(z)XO« d (i?W(z)XOf/X'*i?^(,))^(X'*i?^X')^|, 
1=1 

(3.19) 

where the sup is over all rank two perturbations of X of the form X' = X + xE ki 
where {E k i)ij = 6 ik Sji + SuSij and R' N (z) — (zIn — X'X'*)^ 1 . Additionally, a + 
(3 + 7 + 26 < 4 and each of a, b, c, d, e, /, g, h, q, r are one of i, I, k. The bound (3.14) 
then immediately follows from (3.7) and Lemma 2.4. 

To prove the bound (3.15), we can assume by (iii) of Lemma 2.3 that \x\ < 
eNN- 1 /*. Then 

\\X' N \\ = \\X N \\+o(l). (3.20) 

Additionally, 

R' N (z) = R N (z) + R N (z){xE k iX* N + xX N E tk + x 2 E k iEi k )R' N (z). 

Thus, 

\\R' N (z)\\ < \\R N (z)\\(l + \ X \^\\X N \\) < \\R N (z)\\(l + e N \\X N \\). (3.21) 

Using (3.19-3.21), one obtains (3.15). □ 
It follows from (3.12) and Lemma 3.2 that for i = k 

Ji„ 2 , JWD , , ^/P 4 (|Im(z)|- 1 ) 



(z - a 2 zg N (z) ~ a 2 c N + a 2 )ER u {z) = 1 + O ( ^ ^ n ' ) . (3.22) 



Summing over i and dividing by N gives 

= i+n(-±-l 

N 



( . - a 2 zg N (z) a 2 c N + <j 2 )g N (z) = 1 + ( MjgM^T ) . , , 
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Additionally if i ^ k 

(z - a 2 zg N (z) - a 2 c N + a 2 )E[R tk {z)} = O 



N 



(3.24) 



Now we use (3.22) and (3.23) and the following lemma to complete the proof of 
(3.4). 

Lemma 3.3. On Qn, \gN(z) — a 2 gff(z)\ is uniformly bounded in z and N from 
below by a positive constant. 

Proof. Assume it is not, then for any 5 > there would exist a z such that \qn{z) — 
cr 2 g 2 N {z)\ < 8 2 which in turn implies |ffjv(.z)| < <5 or \a 2 gN(z) — 1| < 6. But 
IsaK z )I < $ contradicts (3.23) if S is sufficiently small once L from (3.8) is chosen 
to make O 



P4(|Im(z)| 
JV 



— ^ small enough because (z — a 2 zg^{z) — a 2 c^ + a 2 ) is 

bounded on Qjv- 

On the other hand if \a 2 gN(z) — 1| < S then 

\(z-a 2 zg N (z)-a 2 c N +a 2 )g N (z)-l\ > \zg N (z)-a 2 zg 2 N (z)\-\o- 2 (l-c N )g N (z)-l\ 

but 

\zg N (z)-a 2 zg 2 N (z)\ < (T + 1)6 

\(a 2 (l - c N )g N (z) - 1) + c N \ = |(1 - c N ){a 2 g N {z) - 1)| < |(1 - c N )5\ 
So for 5 small and L sufficiently large we reach a contradiction with (3.23). □ 

Let s N (z) := 1+<T ' g f (z) 2 ( f then by Lemma 3.3 and (3.22) we have: 



s N (z) - z = 



1 + o- 2 9n(z){c n - 1) 



-z = 



P 4 (|Im(z)|- 1 ) 



(3.25) 



g N (z) - a 2 g%(z) ~ ~ \ N 
Finally, for z € Qn, gn{z) — ga,c N (sN(z)), which can be seen by evaluating 
zo 2 gl. CN {z) + (<j 2 (c n - 1) - z)g a , CN (z) + 1 = 
at Sn(z). This yields: 

1 ^p^^T ^^- Mz)) ~ + a ^ CN ^KMz)) + i = o 

Rearranging this equation and applying the estimates in (3.7) gives: 

2cr 2 a 4 (c N -iy 



(g a .c N (s N (z)) - g N (z)) 



> (9<t,c„(sn(z)) - 9n(z)). 



Ilm(z)!- 1 |Im(z)|- 2 

(3.26) 

So g<T,c N {sN{z)) = gN{z) for sufficiently large z and then on Qjv by analytic con- 
tinuation. 

Combining (3.26) and (3.25) gives: 



cn (^0 dfifj^ CN (x) 
x 



\9N(z)-ga, CN (z) l 

J z - x sn{z) — x 
on Q?q. This completes the proof of (3.4). 



< 



C( SN (z)-z)j^^ 



< 



N 
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Beginning from (3.24) we now finish the proof of (3.5). We have 



g N (z)(z-a zg N {z)-a c + a )E[R ik (z)] =g N (z) + 



P 4 (|Im(z)|- 1 ) 



N 



l + 0{ ^)\-'))) nMik(z)] = /r-M 1 1 



N j J V N 

Recall that z € On, and L can be chosen such that O ^^iiii^MJ — Vj Q n the l.h.s 
of (3.27) is less than 1/2 in absolute value. Then: 

' P 5 (\lm(z)\^) 



E[R lk (z)]=0 



N 



This completes the proof of (3.5). 

Our final step in the proof of Proposition 3.1 is to prove (3.6). First note that if 
zeO r N then 

V[^(z)]< E [||i? N (z)|| 2 ]< i4E[ll ^ (z)l|2] 



\lm(z)\ 4 N 

For the remainder of the proof we will assume z € On- We begin with the resolvent 
identity (3.1) applied to E[R ik (z)Ri k (z)] and then apply the decoupling formula 
(3.2): ' 

N n 

zE[R ik (z)R ik (z)\ = S lk E[R lk (z)]+J2J2 E ^^ R ^ X ^ X ^ ( 3 ' 28 ) 

3 = 1 1=1 

2 

= 5 ik E[R ik {z)} + a 2 E[R ik {z)R ik {z)tv N (X*R N (z)X)\ + ?-E[(R N (z)XX* R N (z)) ik R ik (z)} 



N 

(3.29) 



a 2 



+ —E[(R ik {z)(R N {z)XX*R N (z)) ik + {R N {z)XX*R N {z))iiR kk {z))] + —E[R lk {z)R lk {z)] 

(3.30) 

where tn contains the third cumulant term, p = 2, and the error for truncating at 
p = 2. 

Once again using that Rn(z)XX* = zRn(z) — In gives: 
zE[R lk (z)R lk (z)} = S lk E[R lk (z)} + a 2 zE[R lk (z)R lk (z)tT N (R N (z))] - a 2 E[R lk (z)R tk (z)} 

2 2 

+ ^E[(R N (z)XX*R N (z)) ik Ri k (z)] + ^E[(R ik (z)(R N (z)XX*R N (z)) ik 

+ {R N {z)XX*R N {z)) ii R kk {z))\ + —E[R ik (z)R ik (z)\ + r N 

(3.31) 

Similar to Lemma 3.2 we use the following lemma to complete our variance 
bound. 

Lemma 3.4. For z e C\ M: 

c^w^tr^w)] < mM')\-Wsiz)n (3 . 32) 

< M M^ng [P w (\\X n \\)(\\Rn(z)\\ 2 + \\R N (z)\\^ 2 )] 
rN < ^ (3.33) 
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Proof. The proof follows from the steps taken in the proof of Lemma 3.2. For the 
reader's convenience the third cumulant term is: 

X " r d 2 R t3 (z)X M R lk (z) 



^^ K3 ((A JV ) fei )E 



2iV 3 / 2 

3 = 1 1=1 



1 ™ 

^ /t3 ((^) H )E 



2iV 3 / 2 , 



9 2 (i? w (z)X),, r) d(R N (z)X) il dR ik (z) 8 2 R ik (z) 
tiik{z) + I — — V (Hn(z)A )u- 



(3.34) 

The first subsum is: 
1 ™ 

+ 2(i?(z)X) 4i (i 1 , A r(z)X) fci (^(2)X) fc; i 1 , lfe (z) + 2i? jfe (z)(i? Ar (z)X) fc; i? jfc (z) (3.35) 
+ 2R ik {z)(X*R N {z)) lk {X*R N (z)X) H R ik {z) 
+ 2{R N {z)X)uRkk{z){X*R N {z)X) u Rik{z)) 

The second subsum is: 
1 ™ 

^ K 3 ((4jv)fcl) (2i? lfe (z)(X*i? A r(z)X)„i? jfe (z)(X*i? Ar (z)) ;fe 

+ 2 (P w (i?jv (z)X) kl Ri k (z) (X* R N (z) ) Ifc 

+ 2i? lfe (z)(X*i? 7V (z)X) H (i? A r(z)X) jJ i? fcfc (z) + 2(i? 7 v(z)X) j; (z)(i? 7V (z)X) fe/ (i? A r(z)X) jJ i? fcfc (z) 
+ 2R lk (z)R lk (z)(X*R N (z)) lk + 2R ik {z){R N {z)X) u R kk {z)) 

(3.36) 

The third subsum is: 



2iV 3 / 2 



^72 E«3((^)fei) (2(R N (z)X) u R ik (z)(X*R N (z)X) u R kk (z) 
i=i 

+ 2{R N {z)X)u{R N {z)X)u{R N {z)X) kl R kk {z) 

+ 2(R N (z)X) u R lk (z)R kk (z) + 2{R N {z)X) a R lk (z)(X*R N {-z)) lk {X*R N {z)) lk 
+ 2(R N (z)X) u (R N (z)X) u R kk (z)(X*R N (z)) lk ) 

(3.37) 

Once again by (3.7) and the Cauchy-Schwarz inequality this term is bounded by 
'P 2 (|Im(z)|- 1 )E[P 4 (||X Ar ||)||i? Ar (z)|| 2 ] 



O 



N 



The error term due to the truncation of the decoupling formula at p = 2 is estimated 
as in (3.19-3.21) in Lemma 3.2. □ 

Then using (3.12) to subtract E[Ri k (z)]E[Ri k (z)] from (3.31) gives: 
(z - <j 2 zg N (z) + <7 2 (1- c N ))(E[Ri k (z)Ri k (z)} - E[R lk (z)]E[R lk (z)}) (3.38) 
_ o ^ P 3 (|Im(z)|- 1 )(E[P 10 (||(X jv ||)(||i 1 ' j v(z)|| 2 + p j v(z)|| 3 / 2 ) ^ ^ 
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Repeating the argument for (3.5) leads to: 

-P4(|Im(z)|- 1 )(E[P 10 (||(X Ar ||)(||i? Ar (z)|| 2 + 11^(3)113/2) 



V(Rik(z)) = 



N 



4. Functional Calculus 



□ 



We now extend the results of Section 3 from resolvents to a more general class 
of functions. To do this we use the Helffer-Sjostrand functional calculus ([10], [7]). 
Let / e C' +1 (R), functions with I + 1 continuous derivatives that decay at infinity 
sufficiently fast. Then one can write 

/J *.(«)** • !4(S+'f) <"» 

where: 

i) z = x + iy with i,i/eM; 

ii) f(z) is the extension of the function / defined as follows 

/>ME^«>(v); ( , 2) 

here a e C°°(R) is a nonnegative function equal to 1 for \y\ < 1/2 and 
equal to zero for \y\ > 1. 

From its definition one can see that (4.2) satisfies the following bound: 

y=(x + iy)\ <Constmax(\^j(x)\, 0<j<l + lj \y\ l . (4.3) 

Proposition 4.1. Let be an Nxn real ( complex) matrix that satisfies condition 
C2 (CI). LetM N = ±A N A* N . 

(i) Let f:R^R such that supp(f) nl + C [0, L] for some L > and 
ll/llc 7 ([o.L]) < °°; then there exists a constant, C(L,a,m,4,), such that: 

nf(M N ) u ] - J f(x)d^ CN (x) < C( J L,a,m 4 )||/||c 7 ([o,L])^ 1 for 1 < % < N. 

(4.4) 

(ii) Let f: K — > R such that ||/||c 6 (R + ) < 00 > then there exists a constant, C '(a, mi), 
such that: 

\E[f(M N ) ik }\ < c7(a,™ 4 )||.f||c«(K + )/V- 1 for 1 < i ± k < N. (4.5) 

(Hi) Let f : M — > M such \\f\\ s < oo, for s > 3 £/ien i/iere exists a constant, 
C(s, a, m 4 ) such that: 

\n.f(M N h}\ < C( S , ( 7,m 4 )||/|| s iV- 1 for 1 < i, j < N, (4.6) 

The proof follows as in [17]. We sketch the details below. 

Proof. First, we note that since Mn is a non- negative definite matrix, changing the 
values of f(x) for negative x does not have any effect on the matrix values f(M^)ij. 
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For example, we can always multiply / by a smooth function ip which is identically 
1 on R + and on (—oo, — S\. By the Helffer-Sjostrand functional calculus we have: 



TT J c dz 



dxdy 



df 



- [ ^ (g<r{z) + en) dxdy 

7T J c OZ 



Where 



E[f(M N ) u ] = E 
= E 

= J f{x)dn„,c N (x) — E i J ^eu dxdy 
PedlmWI- 1 ) 



(4.7) 



\E[R u {z) - g a , CN (z)\ < 



N 



by (3.4). Combining this inequality with (4.3), letting I = 6 yields: 



E 



7T J r 8z 



en dxdy 



< qi/llcwi)^" 1 



(4.8) 



Completing the proof of (4.4). The proof of (4.5) follows similarly. 

The rest of the proof of (4.6) follows the proof of Proposition 4.2 in [17], using 
Proposition 1 from [22]. 

We first consider the diagonal case, i = j, without loss of generality let i = 1 
and define the random spectral measure 



N 



H(dx,uj) :=^<5(z-A,)|<Ml)| 2 



(4.9) 



Where A; are the eigenvalues of Mm and <j>i are the corresponding normalized eigen- 
vectors. 

Proposition 2.2 of [17] applied to the measure (4.9) gives 

/•oo />oo 

nf(M N )n}<C s \\f\\ 2 s dye-Vy 2 *- 1 V[R n (x + iy)]dx. 

JO J -oo 

Using (3.6) we can estimate Y[Rij(x + iy)]dx from above by 

P^y- 1 ) " 



N 



-E 



/ 

J — ( 



P 10 (\\X n \\)(\\Rn(x + iy)\\ 2 + \\R N (x + iy)\\ 6 < 2 )dx 



(4.10) 



Once we open the brackets, we obtain two terms. Here, we bound the first term. 
The other term can be estimated in a similar way. 



P w (\\X N \\)\\R N (x + iy)\\ 2 dx 

l|X N A^r|| 



< 



(4.11) 



3 io(||*;v 



y dx + 



\\XnX* 



\x\>\\x N x* N \\ (x- \\X N X* N \\) 2 +y 2 



dx 



< P 2 (t/- 1 )EP 12 (||X Ar | 



(4.12) 
(4.13) 
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By Lemma 2.4 (4.13) can be bounded by CP2(y x ). This leads to 

V[f(X N ) n } < C.Mf J" dye-yy^Peiy- 1 ). 
The integral converges if s > 3. 

In the off-diagonal case i ^ j, we consider the (complex- valued) measure 

N 

/j,(dx,u) :=^2S(x- Xi)4>i(i)(t>i(j), 
i=i 

which is a linear combination of probability measures, and apply Proposition 2.2 
of [17] to each probability measure in the linear combination. Proposition 4.1 is 
proven. □ 



5. Resolvent CLT 

Let m be a fixed positive integer and let i% (z) denote the to x to upper-left 
corner of the resolvent matrix, Rjy(z). Define 

* N (z) = Vn - ga, CN (z)l) ,z e C \ [0,<j 2 (1 + v^) 2 ]- 

Clearly, $n{z) is well defined for z e C \ M. By Lemma 2.2, ^n(z) is well defined 
for z G M\ [0, cr 2 (l + y/c) 2 } with probability going to 1. We are interested in studying 
the random function ^n(z) whose values are in the space of complex symmetric 
to x to matrices. We also define 



<p(z, w) = E 



Vl/c 



w 



Vl/c 



where ■q 1 / c is a Marchenko-Pastur distributed random variable with ratio index i 
and scale index a 2 . We introduce the following notation 



ip ++ (z,w) = E 

f (z, w) = E 

ip^ (z, w) — E 



Re 
Im 

Re 



z - Vi/, 
z 

z - Vi/i 
z 



Re 
■ Im 
• Im 



w 



W - JJi/ c 

w 

W - T] 1/c 

w 



Vl/c 



w 



Vl/c 



Theorem 5.1. Let An be a N x n random matrix with real entries that satisifi.es 
condition C2. Let to be a fixed positive integer and assume that for 1 < i < to 



(5.1) 



exists and for all e > (1.11) holds. Let 

K4,(i) '■— 1714(1) — 3<t 4 , 1 < i < m. 

T/ien the random field ^>n(z) converges in finite- dimensional distributions to the 
random field ^(z) = y/cg 2 c (z)Y(z) where 

Y(z) = (Y ij (z)) 1 ^ m 
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is the Gaussian random field such that the entries Yij(z), i < j and Yu{w), k < I 
are independent when ^ {k,l) and 

Cov(Re Y u (cz),ReY u (cw)) = K 4 (i)Rc[zg aA/c (z)]Re[wg aA/c (w)} + 2a 4 ip ++ i z , w) , 

Cov(lmY u (cz),lmY u (cw)) = K 4 (i)lm[zg aA/c (z)] lm[wg aA / c (w)] + 2a i ip — (z,w) , 

Cov(ReY u (cz),ImY u (cw)) = K 4 (i)Rc[zg aA/c (z)}Im[wg aA/c (w)} + 2a 4 (p+- (z,w) , 

Cov(Re Yij(cz), Re (cw)) = cr 4 (^++ (z, w) ,i 7^ j, 

Cov(lmYij(cz),lmYij(cw)) = a A ^p (z,w) ,i =^ j, 

Cov(Re Yij (cz), lmYij(cw)) = cr 4 ^ (z, w) ,i =/= j. 

In the Hcrmitian case we have the following. 

Theorem 5.2. Let An be a N x n random matrix with complex entries that satisi- 
fies condition CI. Let m be a fixed positive integer and assume that for 1 < i < m 

m 4 (i) := lim -VEU,,! 4 (5.2) 

j 

exists and for all e > (1.11) holds. Let 

n 4 (i) := m 4 (i) — 2er 4 , 1 < i < m. 

Then the random field ^nI?) converges in finite- dimensional distributions to the 
random field ^f(z) = y/cg 2 _ c (z)Y (z) where 

Y(z) = (Y ij (z)) 1 ^ m 

is the Gaussian random field such that the entries Yij(z), i < j and Ym{w), k < I 
are independent when (i,j) 7^ (k, I) and 

Cov(Re Yu(cz),ReYu{cw)) = K 4 (i)Rc[zg (Til/c (z)}Re[wg aA/c {w)} + o"V++ (z,w) , 
Cov(lmY li (cz),lmY li (cw)) = K 4 (i)lm[zg aA/c (z)} lm[wg a ^/ c (w)] + a 4 (p — (z,w) , 
Cov(ReYu(cz),lmY u (cw)) = n 4 (i) Re[zg aA / c (z)] lm[wg a ^/ c (w)] + a i tp+- (z,w) , 

Cov(Re Yij(cz), Re Y i:j {cw)) = ^cr 4 (ip++ (z,w) + ip — (z,«i)),i/j, 
Cov(lmY lJ (cz),lmY ij (cw)) = ^ct 4 (ip ++ (z,w) + ip — {z,w)) ,i ^ j, 

Cov(Re {cz),lmYij{cw)) = ^ct 4 (<£+_ (z,w) - <p+-(w,z)) ,i ^ j. 

Remark 5.3. We remind the reader that the covariance values in Theorems 5.1 and 
5.2 are stated in terms of the Marchcnko-Pastur law with ratio index - and scale 

c 

index a 2 . In some cases it may be more convenient to state the covariances in terms 
of the Marchenko-Pastur law with ratio index c. Indeed, a simple computation 
reveals that for c > and any continuous function /, 

E [f(Vi/c)] = \® [f (^)] + /(°) 

where n c is a Marchenko-Pastur distributed random variable with ratio index c and 
scale index a 2 . In particular, we note that 

9a,l/c{z) = 9aA cz ) + - ^ ( 5 - 3 ) 
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<p(z,w) = -1 



cz 



cw 



C2; — ricCW — r], 

We will need the following lemma for the proof of Theorems 5.1 and 5.2. 
Lemma 5.4. Let B be an N x n matrix. Then 

B*(z - BB*)~ 1 B = B*B(z - B*B)- X 
for all z i Sp(BB*) U {0}. 

Proof. Choose z <£ Sp(BB*) U {0} such that \z\ > \\BB*\\. Then we have that 



(5.4) 



B 



*{z-BB*)~ 1 B = B*- z ^1 + Y^^(BB*) k ^j B 

= B* B \ (l + ^t^B)^ 
= B*B(z — B*B)~ 1 . 

We can now extend the result to all z ^ Sp(BB*) U {0} by analytic continuation of 
the function 

/„„(*) = (B*(z- BB^Bu^v) 
where u, v are arbitrary vectors. □ 

We present the proof of Theorem 5.1 below. The proof in the Hermitian case is 
similar (see also [19] and [17]) and is left to the reader. 



Proof of Theorem 5.1. We write 

/ n \ 



A N = 



4 (m) - 

N — 



\r N J 



( T-m+l \ 

r m +2 



\ r N ) 



where rj is an n- vector representing the i-th row of An- We remind the reader 
that Xn = -kfA^ and we will use the notation X^ — -A=A^™\ Recall that 



Mm = -!rfAMA* M = XnX 



we 



denote by (z) the m x m upper-left corner of the resolvent matrix, Rn(z), of 



Standard linear algebra gives (see e.g. [17] 
R N n) (z)=(zI m -^r i B 



(m) 
JV 



(z)r* 



where 
Let 



B™(z)=I n +X™*(z- 



Y (m) „(m) 



*^ — 1 



X 



(m) 
N ■ 



(5.5) 



r N (z) = 

Then a simple computation reveals that 



^B^\z)r* 



o~ (cn — 1 + zg a ^ CN 
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It will follow from the Central Limit Theorem for Quadratic forms (see the 
appendix of [17]), that ||T jv(^) || is bounded in probability for z e C \ M. Thus, we 
have that 



t>N(z) = gl !CN (z)T N (z) + o(l). 



We note that by (5.5), 



1 



-Tr<>(z) = CAr -l + ^Tr( 



N A W 



(m)*V 
N ) 



and 



^ TYfr y( m ) V-("»)*\-l 



in probability as N — > oo by (3.16). 
Thus, we have that 



= V c w« 

By Lemma 5.4 and (5.5), we have that 
B^\z)=z(z-X^*X^)- 1 = 



(7 

TV 



TW#°(*)J,, 



-TU?™(s)J, 



+ o(l) 
ho(l). 



Z f Z 1 Am)* Am) 

cjv \c N n 



By Theorem 1.4 and (5.6), we note that 



-Tr 



n 
1, 
n 
1, 
n 



ReB^\z)ReB^>(w) 
ImB < ^ l) (z)ImB < ^ l) (w) 
ReB^\z)lmB^\w) 



-Tr 



-Tr 



<P— 
<P+- 



z w 
c c 



/ z w\ 

(--)' 

V c c / 



in probability as TV — > oo. We now claim that for 1 < i < m, 



(5.6) 



-Vk4(4). 3 ReB%?\z) ReB^\w) — > n 4 (i)Rc -g a .i/ c (-) Re —g a ,i/c(—) 

Tl L J 7? L J 77 Lc V C / J L C \ C / 



1 z z 

-Vk 4 (4) 13 ImB^'fz) ImB^'W — ^K 4 (j)Im -ff CT ,i/ c ( - ) 
n ^— ' l j j'i l J l c v c / 



Im 



_ 9<T,1/C , 

c ' v c 



1 Z Z Z/J W 

-Y\ 4 (^iv)ij Re-Bj^z) Im4 m) W — >« 4 (i)Re -&t,i/c (-) Im — 9*,i/c ( — ) 



3=1 
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in probability as N — > oo. Indeed, for the first statement, by the triangle inequality 



Re B { ™\z) ReB^'(w) 



Re 



< 



I Inn 



+ 



Re 



z 

c ' 

z 

-9a,l/. 



Re 



w 

— 9cr,l/c ( 
C ' V C 



(?) 



Imzl 



ReS^H 



-Re 
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W 

— 9a.l/c ( 
C ' V C 



(?) 



By Proposition 3.1, we obtain 



ReB&^z) Re£^' l » 



33 



Re 



-g<7,i/ c , 

c \c 



Re 



+ 0(1). 



The claim is then complete by assumption (5.1). The other two statements follow 
from the same argument. 

Fix p > 1 and consider z\, . . . , z p £ C \ R. We define the family of matrices 



C* 



at 



, 1 < s, t < m 



,(0 



. t and 6g ' t are arbitrary real constants for 1 < s < i < to, 1 < Z < p. We 



where a 

now apply the Central Limit Theorem for Quadratic forms (see the appendix of 
[17]) to the family of matrices C^' and use the above computations to conclude 
that Tn(z) converges in finite dimensional distributions to y^cY(z) for Imz ^ 0. 

For z el\[0,a 2 (l + <Jc)% define 5 = ±dist(z, [0,er 2 (l + y/c) 2 }). Let h(x) be a 
smooth function with compact support where 

h(x) = Qforx(£ [-2<5,er 2 (l + Vc) 2 + 2<5], 

h(x) = 1 forxe [-<5,(T 2 (1 + Vc) 2 + <5]. 

To complete the proof for z G M. \ [0, <r 2 (l + y/c) 2 }, we repeat the same arguments 
as above replacing B^ with ^X^X^B^, where h is a C°°(R) function with 
compact support such that 

h{x) := 1 for x e [-5, cr 2 (l + J~cf + S], (5.7) 

for some 8 > 0. It is essential here that 

F(B% 1) ± h(X N X* N )B^) — ► 

co by Lemma 2.2. 

□ 



as N 



6. Fluctuations of matrix entires for regular functions 
We now prove Theorem 1.6, Theorem 1.7 follows similarly. 
Proof. In Theorem 5.1, Theorem 1.6 is proved for functions of the form 

1 



k 

J2 a i 



i=i 



Zl-X 



-, zi £ [0,er 2 (l + Vc) 2 ], l<l<k. 



(6.1) 



It follows from (1.8) and (5.4) that the limiting variance for functions of the form 
(6.1) given in Theorem 5.1 coincides with the one given in Theorem 1.6. We recall 
that VN(RN(z)ij — 8ijg„ tCN (z)), 1 < i,j < m, converges in finite-dimensional 
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distributions to the random point field v ^ r cg^ c (z)Y i j(z), 1 < i,j < m. For the 
off-diagonal entries i ^ j, one has 

9 lMgl c {w)Cov{Y t3 {z)^{w)) = glAz)gl^w)a\ (£, ^) 

4 2 / x 2 / ^»fc,c(z)-J<r,c(»)) 4 1_c 2 / \ 2 / x 

= fiWl^.cH c ^ _ ^ o- — — 3(TC (z) 5ac H 

= -0- g*A z )9*A w )- [i c - l )9oA w )9oA z ) + — a 2 (z-w) ) ( 6 - 2 ) 

4 1 / s / \ 49a A Z ) - 9a A w ) 

C(J ccr 4 (z — to) 

= ic / 1 1 



C \Z — T] W — T] / 

The calculations in the diagonal case i = j arc similar. To verify that the fourth 
cumulant term in Cov(Yn(z),Yn(w)) gives the required contribution (1.13), one 
uses the identity 



9lA z )9a,i/c (-) = ( z 9lA z ) ~ 9*,c) 

\ C / CO 



1 ({z-a 2 {c-l))g a A z )-l \ 1 {z- ca 2 )gaA z ) - 1 



(6.3) 



Thus, Theorem 1.6 is proved for functions of the form (6.1). By Lemma 2.2, the 
result then also follows for test functions of the form 

k 

Yaihi{x)—^—, Zl $ [0,a 2 (l + ^) 2 }, l<l<k, (6.4) 
tt Zl ~ x 

where hi, 1 < I < k, are C°°(R) functions with compact support that satisfy (5.7). 
By the Stone- Weierstrass theorem (see e.g [21]), functions (6.1) are dense in 

C P (X), for any compact set 1 C K and p = 0, 1, 2, Let / be a C°°(M) function 

such that supp(f) C [— L, L] for some L > 0. Then there exists a sequence of 
functions fj of the form (6.4) that converge to / in the C 4 [— L, L] norm (for the 
definition of the C 4 [— L,L] norm see (1.10)). Then, fj also converges to / in the 
H s norm for s < 4. By the estimate (4.6) 



Y[VN(f(X N ) tJ - f k (X N )ij)] < const\\f(x) - f k {x)\\ s . (6.5) 

Since the r.h.s. can be made arbitrarily small, and oj 2 (f), p 2 (,f) in (1.12), (1.13) are 
continuous in the % s norm, we obtain that y/~N (f(X^)ij — E{f(X^)ij]) converges in 
distribution to the Gaussian random variable defined in Theorem 1.6. Since smooth 
functions with compact support are dense in H s , the result can be extended to an 
arbitrary / G H s . □ 
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